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Abstract 

> . 

The numerical solutions of the non-relativistic Yukawa model on a 3-dimensional size 
. lattice with periodic boundary conditions are obtained. The possibility to extract the 

corresponding - infinite space - low energy parameters and bound state binding energies 
' from eigensates computed at finite lattice size is discussed. 

o 

£ ■ 1 Introduction 

D | The description of deconfined hadronic states from a low energy confining theory is an exciting 
theoretical problem and a challenging numerical task. With the increasing power of computers 
devoted to Lattice QCD and the progress made by algorithms, it is nowadays possible to study 
^ ■ low energy hadron-hadron systems from first QCD principles. 

Contrarily to hadron masses and form factors, scattering properties can not be com- 
puted from infinite volume euclidean simulations pQ . The link between a field theory formulated 
in an euclidean space and the scattering observables is however made possible by studying the 
volume dependence of the bound state spectra in a finite box with periodic boundary condi- 
tions. The underlying formalism was developed by M. Luscher and collaborators in a series of 
papers El HI El El and has been recently updated in view of applications to nuclear and 
hadronic physics H H EDI EI] • 

This formalism was formerly applied to extract the hadron-hadron scattering lengths 
in the quenched approximation [T2], and more recently to obtain the low energy scattering 
parameters in tttt f3J El El El EE] , Kir JH] and even nucleon-nucleon fB] systems from fully 
dynamical lattice QCD calculations. These calculations are crucial in providing some insight on 
the fundamental parameters of the hadron-hadron interaction, as well as low energy observables 
which are poorly known since hardly accessible experimentally. 

*fcsotbor@upo.es 
^carbonel@lpsc. in2p3.fr 
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The aim of our work is to test the applicability of Luscher relations in a model which, 
on one hand, would contain the main physical ingredients of hadron-hadron interaction and, on 
another hand, would be simple enough to be independently controlled. This is provided by the 
Yukawa model, in which a massive scalar particle is exchanged between two fermions. When 
generalized to pseudoscalar and vector exchanges, it constitutes the keystone of baryon-baryon 
interaction models [TU] . 

In the present paper we have considered the non relativistic reduction of the Yukawa 
model that we have solved in a 3-dimensional L 3 lattice with periodic boundary conditions. 
The solutions have been obtained using finite difference schemes in close analogy with the 
methods used in lattice field theory simulations. The L-dependence of the eigenenergies have 
been used to extract the infinite volume low energy parameters, namely the scattering length 
and the effective range. Contrary to the quantum field lattice calculation, these quantities 
can be besides accurately computed by solving the corresponding one-dimensional Schrodinger 
radial equation in an independent way. 

Despite its simplicity, the non relativistic Yukawa model considered here contains an 
essential ingredient of a realistic interaction - its finite range - which plays a relevant role 
in view of extracting the low energy parameters from the finite volume spectra. It offers a 
wieldy and physically sound tool to test the validity of the different approaches discussed in the 
literature, in particular the large and small L-expansion. The full quantum field contents of 
this model have also been considered in a lattice calculation. Preliminary results can be found 
in |23- 

A similar study was undertaken in [§| in the framework of an effective field theory 
without pions EFT (ft) and from a slightly different point of view. These authors fix the values of 
the low energy parameters and obtain the energy levels assuming they satisfy Luscher equations 
[2111312] while in our case all these quantities are independently generated by solving a dynamical 
model and used to find their applicability conditions. 

The properties of the non relativistic Yukawa model in the continuum are briefly intro- 
duced in Section 2. Section 3 is devoted to explain the numerical methods used for its solution 
on a 3-d torus. Results concerning scattering and bound states are displayed in Section 4 and 
some final remarks are given in the conclusions. 



We consider a system of two non relativistic particles interacting by a Yukawa potential of 
strength g and range parameter \x 



2 The model 
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This potential is obtained by Fourier transforming the Born amplitude T# 
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Figure 1: Lowest order two-body amplitude of the Yukawa model (J2J) leading to potential 



where ^ is a Dirac field with mass M, $ a scalar field with mass /i and in which the non 
relativistic reduction - u = 1 q 2 = —q 2 - has been applied [201 |2T]. The same potential results 
from a model with purely scalar fields. 

The reduced S-wave radial Schrodinger equation (in h — c — 1 units) reads 



a 2 e~^ r 

Mo "(r)+M E + ^ 

An r 

and the corresponding tridimensional wavefunction is 

t f. 1 u (r) 



u {r) 







47r r 



(3) 



(4) 



The solutions of © depend a priori on the three parameters (M, /x, g) but written in 
terms of the dimensionless variable x = fir, this equation is equivalent to 



G- 



x 



u(x) = 



(5) 



which depends on a unique dimensionless strength parameter 

g 2 M 



G(M,n,g) 



47T /i 



In particular, the bound state energy (E) and scattering length (a ) of the two-body 
system are given by the scaling relations 



A (G) 



do/1 



(6) 
(7) 



where e(G) and Aq(G) are respectively the binding energy and scattering length of the dimen- 
sionlees problem (J5J). 

We will focus hereafter on finding the solutions of a unit mass particle in the potential 



V(x) = —G- — 



X 



what we will call the dimensionless Yukawa model. All length parameters involved must be 
therefore understood in units of the inverse exchanged mass pT . The lagrangian (J2J) gives 
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rise to an attractive potential, but disregarding its relation with the underling field theory, the 
repulsive case G < can be considered as well. 

The asymptotic norm N s of a bound state with energy e is defined as 

N s = lim u (x)e +k ° x k = (9) 

with Uq normalized by 

POO 

/ ul(x)dx = 1 
Jo 

The convention used for the scattering length is such that the regular zero-energy 
solution of (jSJ) is asymptoticaly given by 



u {x) ~ x + A 

x—>oo 



(10) 



or equivalently the low-energy (S-wave) phase-shifs given by 

S (k) = -A k + o{k 2 ) k 2 = e 
It can be shown that, in the limit of weak coupling constant, one has 

Aq(G) = -G + o(G 2 ) 



which correspond to the Born approximation of the Schrodinger equation with the Yukawa 
model (see Appendix IXjl. 




Figure 2: Bound state binding energies for the ground (e ) and first excited (ei) S-wave state 
as a function of the coupling constant G in the dimensionless Yukawa model (JHJ) 



The "universal" functions e(G) and Aq(G) are displayed in figures 121 and El The system 
has a first bound state (eo) for coupling constant G above some critical value Go ~ 1.680 at 
which Aq(G) has a pole. An infinity of similar branches - and corresponding Aq(G) poles - 
appear for the excited states: e\ at G\ = 6.44, e<i at Gi = 14.34 . . . 
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Figure 3: Scattering length A Q in the dimensionless Yukawa potential (JBJ) as a function of the 
coupling constant G. Dashed line correspond to the repulsive case (G < 0). The tangents 
at the origin indicates the corresponding Born approximation Aq(G) = ±G and the vertical 
asymptote the two first bound state singularities at Go = 1-68 and G\ = 6.44. 



It will be of further interest to consider also the effective range parameter R which 
appears in the low energy expansion 

kcotS (k)=-— + ~R k 2 + o(k i ) (12) 

An 2 

It is worth noticing that the solutions obtained by inserting potential (JH) in the 
Schrodinger equation correspond to summing up the iterations of the diagram displayed in 
figure ^ i.e. the so called "ladder" approximation. This represents a small - though infinite - 
part of the diagrams included in the interaction lagrangian (J2J). Attemps to include all the in- 
teraction terms of this lagrangian in two-body bound and scattering states have been presented 
in [2l] in the framework of lattice quantum field calculation. 

Though the solutions of (jSJ can be easily obtained with great accuracy, we are interested 
on the results provided by the same numerical schemes than the ones used on lattice calculations. 
Thus, by using finite difference methods on an equidistant grid with stepsize a, i.e. 

a 2 it" (x) = u(x + a) — 2u(x) + u(x — a) 

one obtains the ground state binding energies e and low-energy scattering parameters given 
respectively in Tables [T] and 121 

One can remark the different sensibility of the bound and scattering results with respect 
to the stepsize. While low energy parameters vary only of few percent for a G [0.1, 0.5], except 
near the resonant value G = 1.5, the binding energy varies by a more than a factor 2 in the 
same interval. 
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a=0.50 


a=0.20 


a=0.10 


Exact 


G 


eo 


eo 


eo 


eo P 


2.0 


-0.0076 


-0.0177 


-0.0198 


-0.0206 5.63 


2.5 


-0.0791 


-0.1270 


-0.1364 


-0.1397 2.50 


3.0 


-0.2210 


-0.3386 


-0.3625 


-0.3710 1.69 


4.0 


-0.6831 


-1.0665 


-1.1531 


-1.1849 1.07 



Table 1: Binding energy e and radius p = V < x 2 > for the ground state of the Yukawa model 
using finite difference methods with stepsize a. 





a=0.50 


a=0.20 


a=0.10 


Exact 


G 


Aq Rq 


Aq Rq 


A Rq 


Aq Rq 


0.10 
0.20 
0.40 
0.50 
1.00 
1.50 


-0.1029 40.25 
-0.2170 20.24 
-0.4881 10.20 
-0.6518 8.180 
-2.034 4.076 
-8.039 2.628 


-0.1049 40.16 
-0.2227 20.14 
-0.5015 10.10 
-0.6721 8.078 
-2.177 3.958 
-10.81 2.488 


-0.1052 40.15 
-0.2224 20.13 
-0.5035 10.09 
-0.6751 8.063 
-2.199 3.940 
-11.40 2.467 


-0.1053 40.15 
-0.2226 20.12 
-0.5041 10.08 
-0.6761 8.058 
-2.207 3.934 
-11.61 2.460 



Table 2: S-wave low energy scattering parameters (scattering length Aq and effective range Rq 
for different lattice step values a. 



3 Solutions on a torus 

Let us consider now the solutions of the dimensionless Yukawa model (JBJ) on a the 3-d torus, 
i.e. the solutions of 

[-Ag + V L {x)\ $(£) = e(L) $(£) (13) 

satisfying periodic boundary conditions 

$(x + naL) = $(x) \/n E Z 3 

where a > is the lattice step and L e N is the number of lattice points on each spatial 
dimension. 

The "lattice" potential Vl is defined as 

V L (x) = V(x + naL) (14) 

where V is the infinite volume interaction given in (JHJ). Vl incorporates the same periodicity 
than the solutions and contains the interaction with the "surrounding world". 

For large values of aL, Vl is well approximated by V but below aL ~ 10 the "mirror" 
contributions are sizeable and dramatically modify the low-L behaviour of the observables. We 
display in figure 0] the comparison between Vl - converged with respect to the sum over mirror 
images (fTlj) - and V obtained with a = 0.2 and different values of aL. 

In order to obtain the numerical solutions of ()13|) . we introduce in the 3-d cubic lattice 
the coordinates 

= a ( ^ - — j /i = l,2,3 v = 0,1,..., L- 1 
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Figure 4: The lattice potential Vl defined in (114)1 is compared with the infinite volume inter- 
action V for aL = 2 and aL = 10 



Using the finite difference scheme on gets the following expression for the 3-d laplacian 

A$, = -| + (15) 



where % labels the set of coordinate indices i = 12, ^3} and i+fi denotes the nearest neighbour 
of i in the direction \l. 

The stationary states can be obtained by solving directly the eigenvalue equation (fT3*|) 
together with ()15j) . However, in order to follow as close as possible the lattice methods we 
alternatively use an euclidean time-dependent approach. To this aim we consider the evolution 
operator between t and t + A 

U(t,t + A) = e~ iHA 
The time-dependent wave function propagates in this interval according to 

V(x,t + A) = U(t, t + A)*(z, t) = e- iH %e- iH %y(x, t) 

which can be written in the form 

e +i r H y(x,t + A)=e- i % H y(x,t) (16) 

For short A values, we approximate U by its Taylor expansion up to A 2 terms, and get the 
following relation between the wave functions at two consecutifs time steps t and t + 1 

1 iA it 

1 + fH J 1 } 

where we use the notation = \l/(x M (j), tA). 

This approximate numerical scheme has the advantage of preserving the unitarity. 
Using relation ()17|) would however imply the inversion of 1 + iH which is an unpleasant task. 
This can be avoided by noting that ()17|) is actually equivalent to the system of equations 
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-*5 + 4 



The unknown is a solution of an inhomogeneous Schrodinger-like equation which can be 
solved using the same discretization schemes described above for the stationary states. 

The binding energies can be then obtained by propagating an arbitrary initial solution 
rjj in the euclidean time r = it. Indeed, by expanding the initial state in terms of stationary 
eigenfunctions $ ( - n ^ 

n 

it follows that 

W = Y c n $; n) e- e " r A « co$S 0) e- eorA 

n 

where eo is the ground state energy and $^ the corresponding eigenf unction. The stationary 
wavefunction will be normalized according to 

K) 3 E^ 2 = 1 ( 18 ) 




Figure 5: Convergence of the ground state effective energy €q(t) defined in fTT)|) as a function 
of the euclidean time r for two different spatial lattice sizes aL. The values are normalized to 
eo(oo). The ground state energy is well separated from excited states at rA »s aL. Results 
correspond to G — 0.4 a = 0.2 and A = 0.1. 

One can get rid of the spatial degrees of freedom by defining at each time-step r 

This procedure is very similar to the time-slice approach in lattice calculations and can be used 
to disentangle the ground state from the first excitation energies. To this aim it is interesting 



S 



to define the effective energy of the state 

eo(r) - 



1 C(r) 



(19) 



and study the convergence eo( r ) — * e o- The density of states increases with the size of the box, 
making their separation increasingly difficult. However until moderate lattice sizes aL ~ 10 - 
values we will be further interested in - the energy levels are well isolated, what makes easy 
the convergence of the euclidean-time propagation of the ground state. This is illustrated in 
Figure El showing the r-dependence of the ground state effective energy ()19j) - normalized to 
eo(oo) - for two different aL values obtained with G = 0.4, a = 0.2 and A = 0.1. In these cases 
the results are converged for rA rs aL. 

4 Results 
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L 



Figure 6: L-dependence of the first free two-body energy states 



Before presenting the results of the Yukawa model for bound and scattering states we 
will summarize the free case. The non relativistic energies of two - unit mass - particle states 
in a lattice with periodic boundary conditions are given by 



4?W = l±{i 

u=l ^ 



2tt 

— cos — n a 
L M 



4 



y^sin" 



7T 



(20) 



with rifj, — 0,1,2 . . . , L — 1. The continuum limit of an energy state {n^} is reached - indepen- 
dently of a - for L » 7rMax{ri M } with 



2tt^ 
aL 



The ground state energy (n^ = 0) is zero for any value of L and the corresponding 
wavefunction is a constant, which normalized according to (fTH)) . equals 



(0) 



(21) 
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All excited state energies depend on the lattice size and decreases asymptotically like l/(aL) 2 , 
a regime already reached with a number of lattice points L ~ 10 for the first excitation. The 
lowest part of the free spectrum is displayed in figure El One can remark the high degeneracy 
of the states. 




5 10 15 

al_ 

Figure 7: Evolution of the first free e n (L) trajectories (solid lines) with the interaction strength 
in the Yukawa model. Results correspond to G = 0.4 (dotted) and G = 1.5 (dashed), both 
below the bound state threshold and G = 2.5 (dotted-dashed) where the two-body systems has 
a bound state with e = —0.12. 

The generic evolution of the free trajectories when increasing the interaction is dis- 
played in Figure The ground e (L) and first excited €i(L) state energies are plotted as a 
function of the lattice size L for different values of coupling constant G. 

In the attractive case we are considering, the ground state energy eo becomes negative 
and L-dependent. For G < Go, i.e. before the appearance of the first bound state, the excited 
states remain with positive energy and all trajectories tend to zero in the limit (aL) — > oo. 
This behaviour is illustrated for G = 0.4 (dotted line) and G = 1.5 (dashed line) in figure 
In this range all states tend to the free solutions e n (L) — * e£ (L) for large enough values of the 
lattice size. 

When Gq < G < Gi the ground state trajectory eo(L) tends to a non zero negative 
value corresponding to the infinite volume binding energy eo while the first excited state e\(L) 
tends to zero as for the G < G case. Results corresponding to G = 2.5 are displayed in figure 
E|in dot-dashed line. The sign of the e\(L) will in fact depend on the value of G. For G close to 
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the first bound state pole Go (see figure OJ), €i(L) is always positive while for larger values of G, 
close to the second bound state singularity Gi, €i(L) is negative. The change of sign depends 
actually on L but for large enough lattice sizes it corresponds to the zero of A (G). 

It is interesting to remark that despite the singular behaviour of A (G) displayed in 
figure El the G-dependence of e n (L) at a fixed lattice size remains smooth even when crossing 
the bound state singularities. This is illustrated in figure |H1 where the e n (G) dependence is 
shown for two different values of aL. Notice the change of sign of ei(L) at G ~ 4.5 for which 
A (G) « 0. 

A striking difference between the free and the interacting case - not clearly manifested 
in figure 0- is the e n (L) dependence in the limit of small lattice sizes. We will see in the next 
section that for interacting systems, this dependence is ~ ,JjL 3 while for the free case they 

2 

behave like ~ -ttto. 

(aL) 2 

The behaviour described above is qualitatively the same in all the intervals G n < G < 
G n+ i of the coupling constant for the corresponding e n (L) state. 




Figure 8: Dependence of e n (L) on the coupling constant at fixed lattice size. Results correspond 
to aL = 4 (solid line) and aL = 6 (dashed) with a = 0.1. 



Whatever the value of the interaction strength, the two-body spectrum in a finite box 
is always discrete. The main achievement of Luscher and collaborators is to take profit of the 
Independence of e n - or more precisely of the differences e n — ei°^ - to extract the infinite volume 
low energy parameters and bound state energies. Such a possibility will be discussed in the 
next subsections. 
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4.1 Low energy scattering parameters 

For a coupling constant G < G$ = 1.68, the two-body system in the infinite volume has only 
scattering states while in a finite box the spectrum is constituted by a series of discrete values. 

In his first work devoted to this subject [3], Luscher established a relation between 
the two-body binding energy e n (L) on a 3-dimensional spacial lattice with periodic boundary 
conditions and the corresponding scattering length. For the S-wave ground state energy eo(L) 
it reads [22] 




where A is the infinite volume scattering length and C\ = 2.837297 and c<i = 6.375183 are 
universal constants, independent of the details of the particular dynamics. This relation was 
proved to be valid in non relativistic quantum mechanics as well as in quantum field theory and 
must be considered as an asymptotic series on powers of Ao/(aL). For attractive potentials, 
and with our convention for the scattering length, one has Aq(G < Go) < and consequently 
e (L) < 0. 



Equation (|22j) is the most popular of the Luscher relations and has been widely used 
in lattice calculations to extract the value of Aq from a fit to the computed €q(L). We adopt a 
slightly different point of view by constructing from the €q(L) values, quantities tending to Aq 
- the quantity we are interested in - for large values of aL. This merely consists in inverting 

(J22J). 

To this aim it is interesting to consider slowly varying functions and to use - rather 
than e (L) - the combination 

4 0) (L) -L (aL)h (L) (23) 

It tends asymptotically (aL — > oo) to the infinite volume scattering length A and constitutes 
the zero-th order approximation of Luscher expansion which can be written as 

A { °\L)=A |l + Cl (^)+c 2 (^) +...} (24) 

It is possible to get a series A^ (L) of improved values converging towards Aq by solving 
equation (fiHj) truncated at the order n for a fixed value of L. One thus obtains, for instance 




A ( 2) (L) =aLz { ® (25) 
where is a solution of the cubic equation 

2 , 4 0) (£) 



z(c 2 z +c±z + l, 



aL 



This expansion is however of small practical interest for it requires lattice sizes one or 
two order of magnitude larger than the scattering length, as it was already noticed in reference 

Figure |H1 represents the A^ °\l) results for the ground state obtained with G = 0.40 
and a = 0.20. Dashed lines correspond to the solutions with the infinite volume interaction 
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V alone and continuous line with the full interaction Vl in (JHJ). In both cases, A [ q> tends 
indeed asymptotically towards the physical value A = 0.501 (horizontal dotted line) in very 
good agreement with the results given in tabled with the same value of a, but the convergence 
is very slow. Even for computing such a scattering length value, small with respect to the 
lattice sizes, a consequent number of grid points would be required. With aL = 10 only a 
20% accuracy is reached. One can also remark the sizeable effects of the "interactions with the 
surrounding world" below aL pa 10; these contributions - which are absent in the pionless EFT 
considered in jHj - are essential in providing the very smooth variation of A^(L) observed in 
the whole range of aL. 



-0.1 




-10 -5 5 10 

z 



Figure 9: Zero-th order scattering length 
A { q ] = (aL) 3 e /An for G = 0.4 and a = 0.2 
obtained with the infinite volume potential 
V alone (dashed line) and (solid line) with 
the full lattice potential V L (JHJ) 



Figure 10: Wavefunctions corresponding to 
the results of figure |H1 for selected values of L 
and normalized according to (J 18)) 



Corresponding wavefunctions - normalized according to (fPBj) - for L = 40, 60, 100 are 
displayed in Figure El For small values of L they look similar to the bound states ones (see 
next subsection) but their central values depends strongly on the lattice size. Outside the 
interaction region the wave function behaves - according to (fTTH) - like $(a;) ~ 1 + and 
would thus provide an alternative way to extract Aq in case it was numerically accessible. In 
the infinite volume limit they spread over all the lattice with constant amplitude, like in the 
free case. 

The behaviour illustrated in figure El is in fact generic and easy to understand. On one 
hand we have shown (see Appendix [BJ that, under reasonable assumptions, the aL — > limit 
of A^ is finite. For the Yukawa model it gives 

limf (L) = -G (26) 

aL— >0 

On the other hand, the large aL limit is, by construction, 

lim A^(L)=A (G) 

aL— >oo 
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with Aq(G) given in figure H3 For small values of the coupling constant, the L-dependence of 
is almost flat due to the fact that Aq(G) ~ — G but for increasing values of G the presence 
of the bound state pole in figure El makes the variation range increasingly large. 

The inclusion of first and second order corrections in (|24|) - i.e. A^ and A^ - does not 
significantly improve this situation, specially when dealing with large scattering length values. 
This will be illustrated below. On the other hand it is inconvenient to include higher order 
terms since the additional coefficients involved depend on dynamical parameters others than 
A 0. 




Figure 11: Function S(rj) defined by (|2~HJ) 

Equation (|2*2*|) is in fact a large- L expansion of a more general relation established in 
jS] between the two-body phase shifts and the corresponding energy eigenstates on a finite size 
box. For the ground state and using the notations of reference jO] it reads 

k cot S {k) = — — S{rj) (27) 
naL 

where 

(«L) 2 e (L) _ 4 0) 
47T 2 naL 

k 2 = €q(L) and S is a universal function, independent of the interaction model, defined by 

M<a 

SCn) = lim V — 4ttA (28) 

It is represented in figure ITT1 The domain of interest for the ground state is t] < 0. 

It is worth noticing that equation (|2*Tj) is exact for lattice sizes aL > 2R where R is the 
interaction range, i.e. V(x > R) = 0. For interactions of physical interest this regime is reached 
exponentially and independently of the coupling constant and the Aq value. This constitutes a 
remarkable advantage with respect to expansion (|2*2*j) . 
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Using the effective range expansion (fT2*|) one finds the following relation between the 
low energy parameters and the eigenenergies e (L) 

By setting Ro = 0, this relation can be written in the form 

{ 2n\ 2 , ( aL 



eoW = 77 ) S- 1 



-7T- 



\aL J \ A 

where S" -1 denotes the reciprocal function of S. The Luscher expansion (|22|) is recovered by 
developping S~ l for z — > oo in power series of 1/z: S~ 1 (z) = — - + ^ — Z[ pr + . . .. 

Equation (J29~|) suggests the possibility to obtain A as a function of e (L) . . . provided 
Ro is known, which is in general never the case. 

To check the applicability of (J2*§|) in the Yukawa model we fix the Ro value from the 
infinite volume results given on table El and study the convergence of the scattering length thus 
obtained as a function of the lattice size aL. We denote by A (L) the quantity extracted this 
way, i.e.: 

4 S *(i) = ^ ,,.,„ (30) 



§7r R (aL)e (L) — S 



(aL) 2 e a (L) 
TP- 



It generalizes the series of values Aq(L) define above. 

Figure IT2l shows the results for G = 0.40 and a = 0.20, the same parameters than those 
used in figureEl Only results including the fully lattice potential Vl are hereafter shown. Several 
remarks are in order: (i) represents an improvement with respect to A^ for aL > 5, but 
it is significant only in this particular case due to the smallness of A /aL at small values of G 
(ii) Aq S ^ is practically converged at aL ps 5, distance at which remains almost unchanged 
with respect to its value at aL = i.e. to the Born approximation (in) the effective range 
contribution is sizeable: if Rq is neglected - dashed curve obtained setting R = in (j2Hj) - the 
improvement with respect to A^ - and for aL < 5 even to A^ - disappears. 

Similar results are displayed in figure IT31 for a value of the coupling constant G = 1.5 
near the bound state singularity and a lattice space a = 0.10. The infinite volume scattering 
length - with a = 0.10 - is A = —11.4. It corresponds to a = — ~ —16 fm if a pion exchange 
(/i = 0.14 GeV) is supposed to govern the two-body interaction. This value is close to the 
experimental l So nucleon-nucleon scattering length. One can see from this figure that, despite 
the large Aq value, the convergence of A S \L) is the same as in Figure El a few % accuracy is 
already reached at aL w 5. The effective range contribution is here essential: by setting Ro = 
in (jHOj) one would get a result wrong by a factor two (dashed curve) even at aL = 10. Note 
that for small lattices, A °\l) is practically constant - i.e. €q(L) displays a 1/L 3 behaviour - 
but this constant differs from A by one order of magnitude. One can also remark the slow 

(2) 

convergence of A , which, contrary to the G = 0.4 case, provides practically the same results 
than ^ 0) : even at aL = 10 they differ from the asymptotic value Ao by a factor 4. 



The preceding results show the applicability of (j29|) in extracting the scattering length 
Ao for all cases of physical interest, provided that the effective range Ro is taken into account. In 
practice, both values could be determined by fitting the computed finite size energies eo(L) with 
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Figure 12: Different approximations used in calculating the scattering lenght from energy 
eigenstates: defined in (J23J) . obtained by solving the cubic equation (j25J), ^o 5 '' gi yen 
by (|30|) with Rq = (dashed line) and with Rq taken from the infinite volume results on table 
[21 (dotted line). 



the two-parameter function (J29J) for aL > 5. This procedure is however not very comfortable 
for the fitting function f)29j) is defined implicitly. We alternatively propose to determine Aq and 
Rq by computing the energies €q(L) at two different lattice sizes near aL = 5 and solving the 
linear system resulting from (|29|). This gives, e.g. for the effective range 

1 a ( {aL x fe Q {L x ) \ 1 o f (aL 2 ) 2 e (L 2 ) \ 

7i-^ (31) 

7T e (Li) - e (L 2 ) 

The results obtained using this procedure for G — 1.5 are given in table 01 The 
corresponding infinite volume values calculated with a = 0.10 are respectively Aq = —11.4 and 
Ro = 2.47. One sees from these results that a lattice size aL w 6 provides an A value better 
than 3% though the accuracy for Rq is slightly worse. The inclusion of an additional term - 
Pk A - in the low energy epansion (j!2j) does not improve the results. 

It could have some interest to summarize the different results obtained with equation 
(J2HJ) depending on the way Rq is taken into account. This is done in figure ITU for G = 1.5. Solid 
line corresponds to the result of tableEJ- with an averaged value on abscissa aL = 0.5a(Li + L 2 ) 
- where the two parameters Aq and Rq are simultaneously determined. Dotted-dashed and 
dotted lines are the same than in figure H~3l and correspond respectively to Rq taken from Table 
121 and i?o = 0. Notice the different kind of convergence for A^ and the key role of Rq. 
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Figure 13: The same as figure (|12j) for a coupling constant G = 1.5 near the bound state pole. 

(S) 

Aq is practically converged at aL ~ 5, a lattice size smaller than the infinite volume scattering 
length A (G) = 11.4 (horizontal dotted line). Note that for small lattice, A °^ is practically 
constant - i.e. €q(L) displays a 1/L 3 behaviour - but this constant differs from Aq by one order 
of magnitude . 



This method applies well in all the range of coupling constants. We would like to make 
some final comments concerning the extraction of Aq. The usual way of doing so is by first 
fitting the L-dependence of the finite volume energies by a c/(aL) 3 curve and identifying the 
coefficient c with 4ttAq. In our notations, this is equivalent to use the A^ value, which is the 
leading term of Luscher expansion (ffifj) . As we have previously shown, this approach can lead 
to wrong conclusions when using nowadays available lattice sizes, unless the coupling constant 
- and consequently the scattering length - is very small. In this case A is already given by 
the Born approximation. For the Yukawa model - and for small values of the lattice size - it 
actually coincides with — G. 

As one can see from figures [TJandlTJ A °\l) varies very smoothly even for almost 
resonant systems. Thus, when using lattice sizes aL ~ 2 — 3, €q(L) will be well fitted by a c/L 3 
dependence but, as we have shown in these examples, the coefficient c can strongly differ from 
the infinite volume scattering length. Only the use of equation ()29|) at lattice sizes greater than 
the interaction range - aL > 5 in the Yukawa model - could lead to unambiguous extraction 
of the low energy parameters, provided both Aq and Rq are taken into account. 

Finally, we would like to notice that equation (|2*U|) - corresponding to a e (L) ~ ~J0]f 
behaviour of the ground state eigenenergy - depends only on general properties of the interac- 
tion as we have shown in Appendix El This result is totally independent of the Luscher relation 
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aL 2 




A 


4.0 


4.60 


1.9 


-28.8 


4.6 


5.00 


2.1 


-17.2 


5.0 


5.60 


2.3 


-13.9 


5.6 


6.00 


2.4 


-12.5 


6.0 


7.0 


2.5 


-11.7 


7.0 


8.0 


2.6 


-11.3 


8.0 


10.0 


2.6 


-11.3 


10.0 


12.0 


2.6 


-11.4 


Inf. volume 


2.47 


11.4 



Table 3: Low energy parameters extracted from ground state energy at different lattice sizes 
eo(L) by mean of equations (|31|) and (|29|) . They correspond to G = 1.5 and a = 0.1. Last row 
indicates the infinite volume results. 

(}2Tj) and they have even differents domains of applicability, namely aL « 1 and aL » 1. It 
is striking to note that eo(L) has, in both limits, the same L-dependence. 

The results presented above concern the S-wave ground state energy eo(L). For excited 
states, Luscher expansion (j22J) reads 

e n (L) = e^(L) + A n (L) 

with ei 0) (L) given by ()20|) and A n (L) having the same form with differents values for the 
coefficients q. In practice it has the same drawbacks than those described for the ground 
state. Equation (J29|) can be used, in principle, to obtain A from the first excited state energy 
ei(L), now with 77 > 0. At large values of aL, the L-dependence of excited states e n (L) is 
however dominated by the l/(aL) 2 term of free solutions 6n\L), thus making the calculations 
numerically more difficult. It is worth noticing that in the limit aL — > 0, all the excited states 
tend to the same value than eo but this limit is reached at lattice sizes increasingly small (see 
Appendix EI). This behaviour is illustrated in figure ITBI for the case G = 0.4. 

4.2 Bound states 

Bound states are not affected by the "no go" theorem in the euclidean formulation £Q and can 
by directly computed on lattice simulations at finite volume. However, the finite size effects 
are often non negligible in practice and can be controlled by similar expressions than those 
used for scattering states. In addition, the determination of their binding energies in finite size 
lattice calculations is made difficult by the existence of another length scale, independent of 
the interaction range, given by the size of the state. 

In his pioneer work [2] Luscher established the following relation between the ground 
state energy e of a non relativistic two body system and its value e(L) computed in a finite 
size box with periodic boundary conditions j26j : 

e(L) = e - 6iV s 2 t—- (32) 

where k,q = \/—€q is the corresponding wave vector and N s the asymptotic norm, defined in 
(jnj), and determining the large x behaviour of the wavefunction (j3J) 
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Figure 14: Different results obtained with equation (|29|) depending on the value given to 
Rq. Solid line correponds to the results of table |3] with an averaged value on abscissa 
aL = 0.5a(Li + L 2 ) where the two parameters are simultaneously determined. Dotted-dashed 
and dotted correspond respectively to Rq taken from Tableland Rq = 0. 




Figure 15: L-dependence of the ground eo and first energy e± states for G = 0.4 in the limit of 
small lattice sizes. All states converge to the same value, with different slopes depending on 
their quantum number. 



The two-body energy in a box e(L) tends exponentially towards the infinite volume 
bound state binding energy eo- However its decreasing rate kq contains already the required 
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information and can be consequently extracted well before the asymptotic region is reached. 

We first present the results concerning a deeply bound state. They have been obtained 
with G = 4.0 and a = 0.1 and the e(L) dependency is displayed in Figure ITB1 The asymptotic 
value 6q — — 1.16 is in good agreement with the results of tabled with the same lattice spacing 
a = 0.10. Although this value is reached - at 1% accuracy - from aL = 6, the same accuracy 
can be obtained by doing a two-parameter fit with function (J32j) in the region aL e [3, 4], where 
€o(L) is far from being converged. 

o i j 
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Figure 16: L-dependence of the ground state 
binding energy in the dimensionless Yukawa 
model with G = 4 and a = 0.1, Dashed line 
corresponds to one single term V in expan- 
sion ()14|) and solid line to the full interaction 
Vl- Horizontal dotted line indicates the infi- 
nite volume result. 



Figure 17: Wavefunctions corresponding to 
figure HE! for different values of L 



Corresponding wavefunctions, normalized according to (|18|). are displayed in figure IT71 
for different lattice sizes. The periodic boundary conditions give rise to non vanishing values at 
the lattice edges, which are reminiscent of the free ground state solution (}2"Tj) . The characteristic 
length of this state is p = x 2 > = 1.07, smaller than the interaction range R of the potential 
(JBJ). One can thus expect that, like for the scattering states, its properties would be well defined 
at aL ~ 5. At distances greater than the interaction range, the wavefunctions are not sensible 
to the boundary conditions and shows the typical exponential behaviour of the bound states. 
On the contrary at small values of lattice sizes, they look very similar to the scattering states 
displayed in figure El The difference between bound and scattering wavefunctions can be 
formalised by the parameter s(L) = which tends exponentially to (bound) or like 1/ (aL) 
to 1 (scattering) at large aL values 



In analogy with what we have done for scattering states, it is possible to extract both 
eo and N a by writing relation (|H2*jl at two different lattice spacings. The binding energy is thus 
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determined by the solutions of the non linear equation 
Corresponding results for G = 4.0 and a = 0.1 are given in Table |3J 



(33) 



aL\ (1L2 


Kq 


eo 


2.0 3.0 


1.14 


-1.30 


3.0 4.0 


1.09 


-1.19 


4.0 5.0 


1.08 


-1.16 


5.0 6.0 


1.08 


-1.16 


Inf. volume 


1.07 


-1.15 



Table 4: Ground states energy for G = 4.0 and a = 0.1 by solving equations (J33|) at different 
lattice sizes. Infinite volume results are taken from tabled with the same value of a. 



The binding energy and characteristic length of this example roughly correspond to 
the physical case of 4 He nuclei. Indeed, by inserting the e = —1.16 and jjj = 0.15 in equation 
(JHJ) one gets a binding energy in units of nucleon mass E/M = —0.026 while the experimental 
value for A He is Ei He /M = -0.030. 

The situation is less comfortable when dealing with loosely bound states. An example, 
close to the deuteron, is obtained with G = 2.5 for which e = —0.136 and, according to 
E/M = —0.003. Corresponding results are displayed in figures ITH1 and IT91 The convergence of 
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Figure 18: L-dependence of the ground state 
energy for G = 2.5 and a = 0.1. Notations 
are the same than for figure [T01 



Figure 19: Wavefunctions corresponding to 
figure f° r different values of L 



e(L) is reached in this case at aL > 15 due to the larger size of the state, for which p = 2.5. 
It is however possible, by using equation (J3llj) . to reach already a few percent accuracy in the 
range aL G [8, 10]. Results obtained this way are given in Table El 
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5.0 6.0 


-0.131 


1.2 


6.0 8.0 


-0.143 


1.2 


8.0 10.0 


-0.139 


1.2 


10.0 15.0 


-0.137 


1.2 


15.0 20.0 


-0.137 


1.2 


Inf. volume 


-0.136 


1.17 



Table 5: Ground states energy eo and asymptotic norm N s for G = 2.5 and a = 0.1 obtained 
by solving equation at different lattice sizes. 



These examples indicate that the critical lattice size is not given by the interaction 
range but by the size of the bound state itself. It is nevertheless possible by means of equation 
(J22J) to access the binding energies of loosely bound states like deuteron. Notice however, that 
in order to have a few percent accuracy on the binding energy, a consequent number of lattice 
points L ~ 100 would be in this case required. 

The limit aL — > of 6o(L) obtained in Appendix El holds for bound as well as for 
scattering ones: in both cases the eigenenergies display the 1/L 3 law given by equation (|26|). 
The differences in €q(L) appear clearly in the regions where the interaction plays no role, 
which in the Yukawa model is at aL « 5. They are manifested by calculating the quantity 
(aL) 3 eo(L) j (4irG) - which corresponds to A^ /G of the previous section - and comparing the 
results for bound and scattering states. This has been done in figure EH for increasing values 
of G. For the bound state case (solid lines) one can see a cubic divergence corresponding to 
eo term in (J32|) . For scattering states (dashed lines) one observe a linear dependence before 
reaching an horizontal asymptote Aq(G)/G. Both behaviours are easy to disentangle in the 
extreme cases of deeply bound state or small scattering length but it becomes increasingly 
difficult when approaching resonant scattering and/or loosely bound state, unless huge number 
of lattice points are used [2H] . 

5 Conclusion 

We have obtained the solutions of the non relativistic Yukawa model j!} in a 3-dimensional 
lattice with periodic boundary conditions and examined the possibility to extract infinite volume 
low energy parameters - scattering length A and effective range Rq - and bound state binding 
energies from the computed eigenstates e n (L) at finite lattice sizes aL. The eigenergies e n (L), 
and corresponding eigenfunctions, have been calculated using finite difference schemes and 
propagation in the euclidean time, in close analogy with the methods used in lattice field 
theory simulations. The low energy parameters have been independently computed by solving 
the corresponding one-dimensional Schrodinger radial equation. 

We have considered different approximations of the Luscher relations and shown that a 
lattice size of aL 5 - in units of the meson exchanged Compton wave lenght - - is enough to 
determine the low energy parameters at few % accuracy. The method is based in computing the 
energies at two different lattice sizes and solving a linear system of rank two. Such a possibility 
is independent of the scattering length value and applies in the resonant case as well. We have 
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Figure 20: Comparison of the (aL) 3 eo(L) / (AttG) behaviour for bound (solid lines) and scattering 
(dotted lines) states. At aL — > they all tend to —1 as given by (|2l)j) but at large aL, they 
manifest either a cubic divergence or a slow 1/L convergence towards the asymptotic value 
A (G)/G. 

found in particular that the effective range Rq plays a crucial role in determining the scattering 
length value at moderate lattice volumes. 

We have shown that in the limit of small lattice sizes, the L-dependence of eigenen- 
ergies is dominated by an e n (L) = —G/(aL) 3 term, where G is the strength parameter of the 
dimensionles potential (jSJ). This regime is already reached at moderate sizes aL < 2, making 
ambiguous the extraction of the scattering lenght by means of Luscher expansion (|22|) . The 
behaviour is the same for bound and scattering states and applies to a large class of potentials. 

For loosely bound states, the critical lattice size is determined by the spatial extension 
p of the state. We have shown that for the deuteron case it is however possible to get accurate 
values of the binding energy and asymptotic norm with aL = 10. The parameters of deeply 
bound states are well determined, like for the scattering case, with lattice volumes aL m 5. 

The results have been obtained with a non relativistic model, which is justified by the 
small energies involved in the calculations. Despite its simplicity, the model considered contains 
an essential ingredient of the hadron-hadron interaction - its finite range - which plays a relevant 
role in view of extracting the low energy parameters from the finite volume spectra. It offers 
a wieldy and physically sound tool to test the validity of the different approaches discussed in 
the literature to study the low energy scattering of baryon-baryon or meson-baryon systems 
from a lattice simulations in QCD. 

The results presented in this work have been essentially limited to the ground state of 
central attractive interactions, depending only on one paremeter. The method can be easily 
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applied to more involved interactions, like hard core repulsive terms or non central potentials 
leading to coupled channel equations. 
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A Born approximation 

Let us consider the solution of (J3J) in terms of the equivalent Lipmann-Schwinger equation for 
the K-matrix 

K = U + P.V. {UG K} 

where P.V. denotes the principal value. The lowest order in the coupling constant (g 2 ) is given 
by the Born term, i.e. 

K = U 

The K-matrix is related to the S-wave phaseshifts 5q at momentum fco by [2Zj 

1 tan5o(^o) 



K (k ,k ) 



{2it) 2 m R k 



where K denotes the on-sell L = term of the i^-partial wave expansion and mjj the reduced 
mass of the system. In the limit ko — > on gets 

a = (27r) 2 m^ (0,0) 

The lowest order (Born approximation) is thus given by 

a (g) = (2ir) 2 m R U(0, 0) + oQ? 4 ) 

where U must be understood as the k = k' = matrix element of the potential in momentum 
space U(k, k'). For the Yukawa model one has 

U(k, k>) = -r—i q = k- 9 

K ' (2tt) 3 q 2 + fi 2 

and in case of two identical particles with mass M = 2mn 



9 M 

a o{9) = — A 2 +o{g 

Air n z 

In dimensionless units ((7|) it reads 

A (G) — —G + o(G 2 ) 

which is our equation (JTTJ). 
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B Small volume limit 



We would like to study here the zero volume limit of equation (jT3|) with periodic boundary 
conditions and a potential of the form (j!4jl . To this aim, it is useful to develop the periodic 
wave function in a Fourier series: 

and similarly for the potential, 

V L (x) = J2 V * ^ 

By inserting these expresions into equation (|13|) it results the following infinite system of coupled 
linear equations: 



-l) *~ e 



(34) 



n'gZ 3 



The Fourier coefficients of the potential, Vj%, take a rather simple form due to the 
interactions with the " surrounding world" , 



/•aL 

-r^vr / d 3 fe-^V V(x + aLn) 

—y 

(aL) 3 4^ 

n 

—y 



aL paL paL 

dx\ \ dx 2 I dx s e~^ m ' x V(x + aLn) 
Jo Jo 



(aL) 



(aLf 



(ni+l)aL p(n2+l)aL i>(ri3+l)aL 

dyj dy 2 dy 3 e-^™*V($) 



n\aL 



n^aL 



n^aL 



d 3 ye-^™-$V{y) 



(35) 



For the dimensionless Yukawa potential (jHJ) they read: 



4ttG 



4ttG aL 

(aL) 3 l+4 (m«y 



m = (0,0,0) 
rh ^ (0,0,0) 



The leading term in the small volume limit is given by V$. Neglecting other components, the 
system of equations (|3"4"J) decouple into 



and the eigenenergies are given by: 



n 



C n + V s Cf, = , 



(aL) 3 \aL J 



(36) 



They correspond to the continuum limit of the free result (|2()jl shifted by a constant potential 
of depth Vg. 
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One has for the ground state 



e a (L) = (37) 



a result which is equivalent, and prove, equation (J26J). Notice that in the limit aL — > all the 
excited states e n display also the same behaviour and furthermore they all tend to the 

same value 6q(L), independently of n. The lattice size at which they become negative aL = 
are however increasingly small. 
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